National Scale Real-Time Surveillance of SARS-CoV-2 Variants Dynamics by Wastewater Monitoring in Israel.

In this report, we describe a national-scale monitoring of the SARS-CoV-2 (SC-2) variant dynamics in Israel, using multiple-time sampling of 13 wastewater treatment plants. We used a combination of inclusive and selective quantitative PCR assays that specifically identify variants A19/A20 or B.1.1.7 and tested each sample for the presence and relative viral RNA load of each variant. We show that between December 2020 and March 2021, a complete shift in the SC-2 variant circulation was observed, where the B.1.1.7 replaced the A19 in all examined test points. We further show that the normalized viral load (NVL) values and the average new cases per week reached a peak in January 2021 and then decreased gradually in almost all test points, in parallel with the progression of the national vaccination campaign, during February-March 2021. This study demonstrates the importance of monitoring SC-2 variant by using a combination of inclusive and selective PCR tests on a national scale through wastewater sampling, which is far more amendable for high-throughput monitoring compared with sequencing. This approach may be useful for real-time dynamics surveillance of current and future variants, such as the Omicron (BA.1, BA.2) and other variants.

Human metapneumovirus prevalence during 2019-2021 in Israel is influenced by the COVID-19 pandemic.

OBJECTIVES: To compare infection rates and circulating subtypes of human metapneumovirus (hMPV) before (2019-2020) and after the emergence of coronavirus disease 2019 (COVID-19) (2021) in Israel. METHODS: In total, 12,718 respiratory samples were collected from hospitalized patients of all ages during the years 2019 to 2021 at the Sheba Medical Center in Israel and subjected to reverse transcription-polymerase chain reaction analysis. In addition, whole-genome sequencing was performed to characterize the subtypes of hMPV circulating in Israel between 2019 and 2021. RESULTS: A total of 481 samples were found positive for hMPV. Before the emergence of COVID-19, hMPV peaked in winter months and declined thereafter. In sharp contrast, during the COVID-19 pandemic, we observed a delayed peak in hMPV infection cases and higher infection of young children. Viral sequencing showed a shift in the most prevalent circulating hMPV strain from A2b to B1 during the years 2019, 2020, and 2021. CONCLUSION: Compared with the years before the COVID-19 pandemic, in 2021, hMPV mostly affected young children, and the most prevalent circulating subtype shifted from A2b in 2019 to B1.

Specific Detection of SARS-CoV-2 Variants B.1.1.7 (Alpha) and B.1.617.2 (Delta) Using a One-Step Quantitative PCR Assay.

In this report, we describe the development of a reverse transcription-quantitative PCR (RT-qPCR) assay, termed Alpha-Delta assay, which can detect all severe acute respiratory syndrome coronavirus 2 (SC-2) variants and distinguish between the Alpha (B.1.1.7) and Delta (B.1.617.2) variants. The Alpha- and Delta-specific reactions in the assay target mutations that are strongly linked to the target variant. The Alpha reaction targets the D3L substitution in the N gene, and the Delta reaction targets the spike gene 156 to 158 mutations. Additionally, we describe a second Delta-specific assay that we use as a confirmatory test for the Alpha-Delta assay that targets the 119 to 120 deletion in the Orf8 gene. Both reactions have similar sensitivities of 15 to 25 copies per reaction, similar to the sensitivity of commercial SC-2 detection tests. The Alpha-Delta assay and the Orf8119del assay were successfully used to classify clinical samples that were subsequently analyzed by whole-genome sequencing. Lastly, the capability of the Alpha-Delta assay and Orf8119del assay to identify correctly the presence of Delta RNA in wastewater samples was demonstrated. This study provides a rapid, sensitive, and cost-effective tool for detecting and classifying two worldwide dominant SC-2 variants. It also highlights the importance of a timely diagnostic response to the emergence of new SC-2 variants with significant consequences on global health. IMPORTANCE The new assays described herein enable rapid, straightforward, and cost-effective detection of severe acute respiratory syndrome coronavirus 2 (SC-2) with immediate classification of the examined sample as Alpha, Delta, non-Alpha, or non-Delta variant. This is highly important for two main reasons: (i) it provides the scientific and medical community with a novel diagnostic tool to rapidly detect and classify any SC-2 sample of interest as Alpha, Delta, or none and can be applied to both clinical and environmental samples, and (ii) it demonstrates how to respond to the emergence of new variants of concern by developing a variant-specific assay. Such assays should improve our preparedness and adjust the diagnostic capacity to serve clinical, epidemiological, and research needs.

HiSpike Method for High-Throughput Cost Effective Sequencing of the SARS-CoV-2 Spike Gene

The changing nature of the SARS-CoV-2 pandemic poses unprecedented challenges to the world's health systems. Emerging spike gene variants jeopardize global efforts to produce immunity and reduce morbidity and mortality. These challenges require effective real-time genomic surveillance solutions that the medical community can quickly adopt. The SARS-CoV-2 spike protein mediates host receptor recognition and entry into the cell and is susceptible to generation of variants with increased transmissibility and pathogenicity. The spike protein is the primary target of neutralizing antibodies in COVID-19 patients and the most common antigen for induction of effective vaccine immunity. Tight monitoring of spike protein gene variants is key to mitigating COVID-19 spread and generation of vaccine escape mutants. Currently, SARS-CoV-2 sequencing methods are labor intensive and expensive. When sequence demands are high sequencing resources are quickly exhausted. Consequently, most SARS-CoV-2 strains are sequenced in only a few developed countries and rarely in developing regions. This poses the risk that undetected, dangerous variants will emerge. In this work, we present HiSpike, a method for high-throughput cost effective targeted next generation sequencing of the spike gene. This simple three-step method can be completed in < 30 h, can sequence 10-fold more samples compared to conventional methods and at a fraction of their cost. HiSpike has been validated in Israel, and has identified multiple spike variants from real-time field samples including Alpha, Beta, Delta and the emerging Omicron variants. HiSpike provides affordable sequencing options to help laboratories conserve resources for widespread high-throughput, near real-time monitoring of spike gene variants.

The SARS-CoV-2 Lambda variant and its neutralisation efficiency following vaccination with Comirnaty, Israel, April to June 2021.

The SARS-CoV-2 Lambda (Pango lineage designation C.37) variant of interest, initially identified in Peru, has spread to additional countries. First detected in Israel in April 2021 following importations from Argentina and several European countries, the Lambda variant infected 18 individuals belonging to two main transmission chains without further spread. Micro-neutralisation assays following Comirnaty (BNT162b2 mRNA, BioNTech-Pfizer) vaccination demonstrated a significant 1.6-fold reduction in neutralising titres compared with the wild type virus, suggesting increased susceptibility of vaccinated individuals to infection.

Genomic variation and epidemiology of SARS-CoV-2 importation and early circulation in Israel.

Severe acute respiratory disease coronavirus 2 (SARS-CoV-2) which causes corona virus disease (COVID-19) was first identified in Wuhan, China in December 2019 and has since led to a global pandemic. Importations of SARS-CoV-2 to Israel in late February from multiple countries initiated a rapid outbreak across the country. In this study, SARS-CoV-2 whole genomes were sequenced from 59 imported samples with a recorded country of importation and 101 early circulating samples in February to mid-March 2020 and analyzed to infer clades and mutational patterns with additional sequences identified Israel available in public databases. Recorded importations in February to mid-March, mostly from Europe, led to multiple transmissions in all districts in Israel. Although all SARS-CoV-2 defined clades were imported, clade 20C became the dominating clade in the circulating samples. Identification of novel, frequently altered mutated positions correlating with clade-defining positions provide data for surveillance of this evolving pandemic and spread of specific clades of this virus. SARS-CoV-2 continues to spread and mutate in Israel and across the globe. With economy and travel resuming, surveillance of clades and accumulating mutations is crucial for understanding its evolution and spread patterns and may aid in decision making concerning public health issues.

Rapid and High-Throughput Reverse Transcriptase Quantitative PCR (RT-qPCR) Assay for Identification and Differentiation between SARS-CoV-2 Variants B.1.1.7 and B.1.351.

Emerging SARS-CoV-2 (SC-2) variants with increased infectivity and vaccine resistance are of major concern. Rapid identification of such variants is important for the public health decision making and to provide valuable data for epidemiological and policy decision making. We developed a multiplex reverse transcriptase quantitative PCR (RT-qPCR) assay that can specifically identify and differentiate between the emerging B.1.1.7 and B.1.351 SC-2 variants. In a single assay, we combined four reactions-one that detects SC-2 RNA independently of the strain, one that detects the D3L mutation, which is specific to variant B.1.1.7, one that detects the 242 to 244 deletion, which is specific to variant B.1.351, and the fourth reaction, which identifies the human RNAseP gene, serving as an endogenous control for RNA extraction integrity. We show that the strain-specific reactions target mutations that are strongly associated with the target variants and not with other major known variants. The assay's specificity was tested against a panel of respiratory pathogens (n = 16), showing high specificity toward SC-2 RNA. The assay's sensitivity was assessed using both in vitro transcribed RNA and clinical samples and was determined to be between 20 and 40 viral RNA copies per reaction. The assay performance was corroborated with Sanger and whole-genome sequencing, showing complete agreement with the sequencing results. The new assay is currently implemented in the routine diagnostic work at the Central Virology Laboratory, and may be used in other laboratories to facilitate the diagnosis of these major worldwide-circulating SC-2 variants. IMPORTANCE This study describes the design and utilization of a multiplex reverse transcriptase quantitative PCR (RT-qPCR) to identify SARS-COV-2 (SC2) RNA in general and, specifically, to detect whether it is of lineage B.1.1.7 or B.1.351. Implementation of this method in diagnostic and research laboratories worldwide may help the efforts to contain the COVID-19 pandemic. The method can be easily scaled up and be used in high-throughput laboratories, as well as small ones. In addition to immediate help in diagnostic efforts, this method may also help in epidemiological studies focused on the spread of emerging SC-2 lineages.

The Rise and Fall of a Local SARS-CoV-2 Variant with the Spike Protein Mutation L452R.

Emerging SARS-CoV-2 variants may threaten global vaccination efforts and the awaited reduction in outbreak burden. In this study, we report a novel variant carrying the L452R mutation that emerged from a local B.1.362 lineage, B.1.362+L452R. The L452R mutation is associated with the Delta and Epsilon variants and was shown to cause increased infection and reduction in neutralization in pseudoviruses. Indeed, the B.1.362+L452R variant demonstrated a X4-fold reduction in neutralization capacity of sera from BNT162b2-vaccinated individuals compared to a wild-type strain. The variant infected 270 individuals in Israel between December 2020 and March 2021, until diminishing due to the gain in dominance of the Alpha variant in February 2021. This study demonstrates an independent, local emergence of a variant carrying a critical mutation, L452R, which may have the potential of becoming a variant of concern and emphasizes the importance of routine surveillance and detection of novel variants among efforts undertaken to prevent further disease spread.

Prolonged detection of complete viral genomes demonstrated by SARS-CoV-2 sequencing of serial respiratory specimens.

Accurate and timely diagnosis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is clinically essential, and is required also to monitor confirmed cases aiming to prevent further spread. Positive real-time PCR results at late time points following initial diagnosis may be clinically misleading as this methodology cannot account for the infection capabilities and the existence of whole genome sequences. In this study, 47 serial respiratory samples were tested by Allplex-nCoV test (Seegene), a triplex of three assays targeting the SARS-CoV-2 RdRP, E and N genes and subsequently assessed by next generation sequencing (NGS). COVID19 patients were tested at an early stage of the disease, when all these viral gene targets were positive, and at an advanced stage, when only the N gene target was positive in the Allplex-nCoV test. The corresponding NGS results showed the presence of complete viral genome copies at both early and advanced stages of the disease, although the total number of mapped sequences was lower in samples from advanced disease stages. We conclude that reduced viral transmission at this late disease stage may result from the low quantities of complete viral sequences and not solely from transcription favoring the N gene.

A Unique SARS-CoV-2 Spike Protein P681H Variant Detected in Israel.

The routine detection, surveillance, and reporting of novel SARS-CoV-2 variants is crucial, as these threaten to hinder global vaccination efforts. Herein we report a novel local variant with a non-synonymous mutation in the spike (S) protein P681H. This local Israeli variant was not associated with a higher infection rate or higher prevalence. Furthermore, the local variant was successfully neutralized by sera from fully vaccinated individuals at a comparable level to the B.1.1.7 variant and an Israel wild-type strain. While it is not a variant of concern, routine monitoring by sequencing is still required.

Detection of SARS-CoV-2 variants by genomic analysis of wastewater samples in Israel.

Investigation of SARS-CoV-2 spread and identification of variants in sewers has been demonstrated to accurately detect prevalence of viral strains and is advantageous to clinical sampling in population catchment size. Herein, we utilized an established nationwide system of wastewater sampling and viral concentration approaches to perform large-scale surveillance of SARS-CoV-2 variants in nine different locations across Israel that were sampled from August 2020 to February 2021 and sequenced (n = 58). Viral sequences obtained from the wastewater samples had high coverages of the genome, and mutation analyses successfully identified the penetration of the B.1.1.7 variant into Israel in December 2020 in the central and north regions, and its spread into additional regions in January and February 2021, corresponding with clinical sampling results. Moreover, the wastewater analysis identified the B.1.1.7 variant in December 2020 in regions in which non-sufficient clinical sampling was available. Other variants of concern examined, including P.1 (Brazil/Manaus), B.1.429 (USA/California), B.1.526 (USA/New York), A.23.1 (Uganda) and B.1.525 (Unknown origin), did not show consistently elevated frequencies. This study exemplifies that surveillance by sewage is a robust approach which allows to monitor the diversity of SARS-CoV-2 strains circulating in the community. Most importantly, this approach can pre-identify the emergence of epidemiologically or clinically relevant mutations/variants, aiding in public health decision making.

Comprehensive Analyses of SARS-CoV-2 Transmission in a Public Health Virology Laboratory.

SARS-CoV-2 has become a major global concern as of December 2019, particularly affecting healthcare workers. As person-to-person transmission is airborne, crowded closed spaces have high potential for rapid virus spread, especially early in the pandemic when social distancing and mask wearing were not mandatory. This retrospective study thoroughly investigates a small-scale SARS-CoV-2 outbreak in Israel's central virology laboratory (ICVL) in mid-March 2020, in which six staff members and two related family members were infected. Suspicions regarding infection by contaminated surfaces in ICVL facilities were nullified by SARS-CoV-2 negative real time polymerase chain reaction (PCR) of work surfaces swipe tests. Complete SARS-CoV-2 genomes were sequenced and mutation analyses showed inclusion of all samples to clades 20B and 20C, possessing the spike mutation D614G. Phylogenetic analysis clarified transmission events, confirming S1 as having infected at least three other staff members and refuting the association of a staff member's infected spouse with the ICVL transmission cluster. Finally, serology tests exhibited IgG and IgA antibodies in all infected individuals and revealed the occurrence of asymptomatic infections in additional staff members. This study demonstrates the advantages of molecular epidemiology in elucidating transmission events and exemplifies the importance of good laboratory practice, distancing and mask wearing in preventing SARS-CoV-2 spread, specifically in healthcare facilities.